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ABSTRACT 

We report on a new Chandra exposure of PSR Jl 809— 2332, the recently discovered pulsar powering the 
bright EGRET source 3EG Jl 809— 2328. By registration of field X-ray sources in an archival exposure, we 
measure a significant proper motion for the pulsar point source over an « 11 year baseline. The shift of 
0.30 ± 0.06" (at PA = 153.3 ± 18.4) supports an association with proposed SNR parent G7.5-1.7. Spectral 
analysis of diffuse emission in the region also supports the interpretation as a hard wind nebula trail pointing 
back toward the SNR. 

Subject headings: pulsars: individual (PSR Jl 809— 2332) - X-rays: general 



1. INTRODUCTION 

3EG Jl 809-2328 was one of the brightest unidenti- 
fied hard- spec trum gamma-ray sources detected by EGRET 
(lHartmanetal.1 119991) . ASCA observations o f the EGRET 
error box revealed an extended X-ray source (Robert set al.l 
120011) . suggestive of a pulsar wind nebula (PWN). A 
9.7 ks Chandra ACIS exposure in August 2000 dBraje et al.1 
2002) revealed a point source connected to the non- 
ther mal X-ray/radio nebula, bo lstering the PWN identifica- 
tion. [Roberts & Broganl (120081) then described G7.5-1.7, a 
partial ~ 0.5 — 0.8° radius radio shell and possible super- 
nova remnant; the center of this SNR candidate lies 0.3° from 
the point source, in the general direction defined by the PWN 
trail. 

The PSR/PWN nature of the source was confirmed by the 
discovery of PSR J1809-2332, a P = 147ms, r c = 68 ky 
pulsar discovered by blind sear ch of the Fermi L arge Area 
Telescope (LAT) 7-ray photons dAbdo et al.ll2009h . This en- 
ergetic = 4.3 x 10 35 erg s _1 ) pulsar powers the 7-ray and 
PWN emission and has a timing pos ition consisten t with the 
X-ray point source at the PWN apex (Ray et al]|2011l) . As for 
most 7-ray selected pulsars, we l ack a dispersion m easure es- 
timate of the distance. However JOka et a l. (1999) suggested 
that the diffuse X-ray emission is anti-correlated with molec- 
ular gas in the Lynds 227 dark cloud; if associated this implies 
a plausible d ~ 1.8 kpc. 

To test the SNR association and to improve spectral and 
morphological measurements of the diffuse emission, we 
have obtained a new Chandra ACIS exposure. By carefully 
matching to the original pointing, we have minimized sys- 
tematic effects and allowed excellent frame-referencing and 
astrometry of the point source. The ~ 11 y baseline between 
the exposures thus allows a sensitive study of the proper mo- 
tion and origin of PSR Jl 809-2332. 

2. X-RAY DATA ANALYSIS 

PSR J1809-2332 was observed on July 28, 2011 (obsID 
12546) for 29.7 ks in timed VFAINT mode. As for our Au- 
gust 16, 2000 exposure (obsID 739) we used the ACIS-I array. 
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To minimize changes in the reference star PSFs (see below) 
we matched the roll angle to the 2000 epoch. The point source 
position was not, of course, known in the initial exposure and 
fell rather close to a chip boundary. We thus shifted the point- 
ing ~ 30" in each coordinate to move the point source further 
onto the 13 chip and improve measurement of diffuse emis- 
sion in its vicinity. The data were analyzed with CI AO 4.3, 
after re-processing both frames to the most current ACIS cal- 
ibration. No flaring occurred during obsID 12546, so the full 
29.7 ks was used. For obsID 739, the full 9.7 ks was also us- 
able. 

After combining, exposure correcting (with fluximage) and 
adaptive smoothing, the ACIS-I frame gives a good overall 
view of the pulsar vicinity (Figure 1). There is appreciable 
diffuse X-ray emission across the center of the image. The 
pulsar point source lies near the apex of a patch of brighter 
nebulosity that trails off to the North- West for ~ 5'. On 
smaller scales the PWN has a bright core behind the pul- 
sar, ~ 3 — 5" extension transverse to this trail at the point 
source and fainter diffuse emission 'ahead' of the pulsar. To 
the South of the pulsar, numerous bright X-ray sources from 
the young star cluster S32 are seen. These bright stars provide 
excellent sources to aid frame registration between epochs. 
There are some low significance X- ray co ncentrations outside 
of the bright inner nebula in Figure \Ub)\ None have obvious 
stellar counterparts so it is unclear if these are background 
sources or clumps in the diffuse nebular emission. 

3. FRAME POINTING OFFSET 

To best constrain any pulsar proper motion, we must opti- 
mally register the two observation epochs, using field point 
sources. We identify these using the CI AO wavdetect tool, 
selecting sources with a significance >5 in both observations. 
The young S32 star cluster ~ 8' due South of the pulsar pro- 
vides the bulk of the reference stars; inevitably several bright 
stars fall into chip gaps in one or the other of the observations, 
precluding their use. This leaves 14 stellar sources, plus the 
pulsar. We label the reference stars by increasing phy sical 
pixel x-coordinate (decreasing RA), mark these in Figure [T(a)| 
and tabulate the wavdetect pixel positions and significances 
for each epoch in Table [T] 

The wavdetect tool uses "Mexican hat" wavelet functions 
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Fig. 1 . — Left: Merged Chandra exposure corrected 0.5 — 7 keV ACIS-I array image of the extended emission around PSR J1809— 2332 (observation ID 739, 
12546), adaptively smoothed (0.5 to 15 pixel Gaussian kernels, 10 _6 erg cm _2 s _1 flux minimum). A possible bowshock and trail are seen extending northwest 
from the pulsar. Point sources used in frame matching are shown in green, with the pulsar position marked with a black cross. The white ellipse denotes the 
lcr error ellipse for the extrapolated pulsar position some 10 kyr in the past, while the white arrow points to the proposed explosion site for SNR G7.5— 1.7. 
These two directions agree within la. Right: Merged adaptively smoothed (0.5 to 15 pixel Gaussian kernels, 3 count minimum) 0.5-7 keV image of the 100" 
surrounding the pulsar. Both panels are displayed with a logarithmic stretch. 



TABLE 1 
Wavdetect Source Positions 



Src ObslX ObslY Obsl Sig Obs2 X Obs2 Y Obs2 Sig 



1 


3727.7 ±0.2 


3830.9 ±0.3 


8.9 


2 


3894.4 ±0.2 


3527.3 ±0.3 


35.4 


3 


3971.2 ±0.7 


3213.6 ±0.8 


5.9 


4 


4026.8 ±0.4 


4628.9 ±0.5 


5.4 


5 


4035.3 ±0.4 


3225.6 ±0.4 


24.4 


6 


4049.1 ±0.5 


3177.4 ±0.4 


24.1 


7 


4068.0 ±0.2 


3691.5 ±0.3 


6.1 


8 


4140.3 ±0.2 


3597.1 ±0.2 


31.7 


9 


4144.9 ±0.9 


3194.2 ±0.9 


6.0 


10 


4149.2 ±0.5 


3125.5 ±0.5 


17.8 


11 


4153.8 ±1.0 


3244.0 ± 0.9 


5.7 


12 


4176.4 ±1.1 


3327.8 ± 1.0 


5.2 


13 


4289.4 ± 1.2 


3108.4 ±1.0 


7.3 


14 


4922.8 ±0.9 


4702.8 ±1.1 


7.6 


psr 


4179.65 ±0.12 


4143.61 ±0.13 


29.5 



of different scales to generate a source list and estimate po- 
sition. This approximation may be inadequate far from the 
ACIS aimpoint where the mirror PSF induces substantial and 
systematic changes in the point source count distribution. By 
approximately matching the pointing and roll angle we have 
kept the change in such distortions minimal between the two 
frames. Nevertheless the PSF shape changes across the frame 
are very large and the majority of our reference stars lie > 6' 
from the aimpoint and pulsar. We have improved the astro- 
metric accuracy and further reduced the effect of PSF vari- 
ation by creating a simulated PSF for each point source and 
using these to fit the source position. 

We extract the spectrum of each source with the specextract 
function, fitting the stellar point sources to an absorbed Mekal 



3796.6 ±0.1 


3892.5 ±0.1 


34.3 


3963.4 ±0.3 


3587.3 ±0.4 


6.0 


4044.7 ±0.5 


3272.7 ±0.6 


11.6 


4095.9 ±0.4 


4689.4 ± 0.3 


11.9 


4104.2 ±0.2 


3286.5 ±0.2 


42.4 


4116.4 ±0.4 


3238.1 ±0.5 


15.7 


4137.6 ±0.2 


3752.1 ±0.2 


11.8 


4209.4 ±0.1 


3657.8 ±0.1 


51.0 


4215.1 ±0.5 


3254.8 ±0.8 


6.0 


4216.1 ±0.7 


3187.9 ±0.8 


8.2 


4225.9 ±0.6 


3303.8 ±0.8 


6.1 


4246.8 ±0.6 


3389.5 ±0.4 


8.5 


4360.0 ±0.5 


3168.4 ±0.9 


5.2 


4991.8 ±0.7 


4762.0 ±0.5 


17.7 


4248.44 ± 0.06 


4204.07 ± 0.06 


59.1 



plasma model, and save this spectrum within the Sherpa plot- 
ting program with the save— chart— spectrum command. Af- 
ter re-scaling by 100 x to improve PSF modeling statistics, 
the spectra are passed to CHaRT, which is a web interface to 
the SAOsac raytrace code. CHaRT also takes as input the ex- 
posure time as well as the source position, which is assumed 
to be the wavdetect position. The output of CHaRT is sub- 
sequently input into MARX to create model PSF events files 
at each source position, corrected for the offset of the science 
instrument module (SIM) from the nominal location. 

The modeled source events are binned to 1/8 of the na- 
tive ACIS pixel resolution and recorded in 320x320 image 
files. These are compared with data cut-outs with matched 

2 http ://cxc . harvard . edu/ chart/runchart . html 
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(a) 



(b) 



(c) 



Fig. 2. — (a) Chandra events cutout of Source 8 (2011 epoch), (b) Simulated PSF at this source position, (c) Computed FoM map zoomed 10 x, the green 
(5FoM = 0.5 contour is shown. The extrema in x and y of this ellipse give the positional uncertainty in the source position. 



40x40 native-pixel images. PSF models and data cutouts are 
prepared for each star in both frames. We use these to solve 
for PSF vs. data position shifts. The fitting program uses a 
maximum likelihood 'figure of merit' (FoM, the negative In 
Poisson probability of getting the observed data counts from 
the model PSF), where we run a grid of trial positions shifted 
at the 1/8 -pixel PSF model resolutions and sum the model 
counts to the observed ACIS pixel. The result is a map of 
source likelihood in PSF position shift x and y (relative to the 
nominal wavdetect location). An example of the cut-out, PSF 
model and shift map are shown in Figure 2. 

In order to compute the best-fit position of the source, we 
fit an elliptic paraboloid to the FoM surface. The minimum of 
this surface determines the best fit offset, and the offset error 
is estimated from the £FoM contours. As with Source 8, sev- 
eral of the ellipses had large axis ratios. However an attempt 
to improve registration by fitting along major and minor axis 
did not reduce the scatter. Thus to be conservative, we project 
the ellipses to the x and y pixel axes and fit in this unrotated 
space. We determine that la (wavdetect) errors correspond 
to a SFoM = 0.5, and thus ascribe la error ellipses to the 
region enclosed by this FoM increase above the fit minimum. 
Comparison with sources near the aimpoint, where a circular 
Mexican hat should be an adequate approximation, confirm 
that this provides a good error estimate. A typical fit and un- 
certainty are shown by the green ellipse in Figure 2. 

Table 2 gives the PSF-fit computed pixel offset of each 
source from the wavdetect position, along with estimates for 
the position errors. The mean offset of all reference sources 
(« 0.05 pixels) is not significantly different from zero in ei- 
ther axis, showing that the wavdetect solutions do not impose 
any large systematic offset between our matched exposures. 

We can now use our stellar position estimates to register the 
frames. For the wavdetect positions we use the values in Table 
Q] directly. For PSF-fits we add the additional offsets of Table 
2 and use PSF localization errors. For each source i, we de- 
fine, e.g., Xi as the difference in x coordinate between old and 
new exposure; epoch positional errors are added in quadrature 
to give 5xi . The best fit frame shift is determined by minimiz- 
ing the x 2 — ^Zitii( x i ~ S x )/Sxi] 2 between frames, where 
S x is the frame shift. Similar values are computed for the y 
coordinate. 

Inspection of the tables shows that source 3 has both a large 



TABLE 2 

PSF Fit Offsets from Wavdetect Positions 



Src 



Obsl (X, Y) 



Obs2 (X, Y) 



1 


(-0. 


,17 


±0. 


,26, 


+0, 


.53 ±0. 


.38) 


(+0. 


,08 ±0. 


,12, 


-0, 


,18 


±0, 


.15) 


2 


(-0. 


,01 


±0. 


,18, 


-0. 


.16 ±0. 


.24) 


(-0. 


,22 ±0. 


,52, 


+0. 


,28 


±0, 


.67) 


3 


(+!■ 


,61 


±1, 


,68, 


-1. 


.59 ±1. 


.46) 


(-0. 


,09 ±0. 


44, 


-2 


,05 


±0, 


.85) 


4 


(+0, 


,22 


±0. 


,70, 


+0. 


.19 ±0. 


.62) 


(+0. 


.30 ±0. 


.46, 


+0. 


,35 


±0, 


.37) 


5 


(+0, 


,03 


±0. 


,24, 


-0. 


,26 ±0. 


.33) 


(+0. 


,38 ±0. 


,18, 


-0, 


,09 


±0, 


.22) 


6 


(-0. 


,61 


±0, 


,37, 


-0, 


,24 ±0. 


.39) 


(+0. 


,62 ±0. 


.51, 


-0, 


,24 


±0, 


.49) 


7 


(+0. 


,50 


±0. 


,35, 


-0. 


,77 ±0. 


.44) 


(-0. 


,34 ±0. 


,28, 


-0, 


,13 


±0. 


.32) 


8 


(-0. 


,09 


±0. 


,20, 


-0. 


.06 ±0. 


.23) 


(+0. 


.05 ±0. 


,13, 


+0. 


,11 


±0, 


.14) 


9 


(+!■ 


,36 


±0. 


,24, 


+0. 


,07 ±0 


.36) 


(-0. 


.39 ±0. 


.29, 


+0. 


,17 


±0, 


.33) 


10 


(-0. 


,55 


±0. 


,36, 


+0. 


.10 ±0. 


.46) 


(+1. 


,83 ±0. 


,85, 


-1 


,46 


±0, 


.83) 


11 


(+0. 


,65 


±0. 


,33, 


-0. 


.51 ±0. 


.54) 


(-0. 


.96 ±0. 


,32, 


+1. 


,01 


±0, 


.28) 


12 


(-0. 


,38 


±0, 


,73, 


+0, 


,68 ±0 


.82) 


(-0. 


,96 ±0. 


.54, 


-0, 


,46 


±0, 


.45) 


13 


(-1. 


,20 


±0, 


,92, 


+1. 


,40 ±0 


.96) 


(-0. 


,98 ±0. 


,64, 


+1. 


,24 


±0, 


.60) 


14 


(+!■ 


,16 


±0. 


,37, 


-1. 


.05 ±0 


.39) 


(+0. 


,98 ±0. 


,27, 


+0. 


,02 


±0, 


.21) 


psr 


(-0. 


,16 


±0. 


,12, 


-0. 


,11 ±0 


.12) 


(+0. 


,04 ±0. 


.06, 


-0, 


,15 


±0, 


,06) 



shift between frames and a large PSF-fit error. We suspect 
that it may be a confused double with variable components 
(hence the poor PSF fit and large apparent shift). We thus 
omit this star from the registration, leaving 13 field stars. We 
also attempted recursive pruning of field stars with the largest 
remaining Sxi ja Xi or 8yi / a Vi ; as expected the nominal fit er- 
rors decreased slightly (especially for the PSF-fit solution). 
However the shifts were small and decreasing the number of 
fit stars may increase sensitivity to systematics. We thus con- 
servatively retain all stars except the obvious outlier source 3 
and report the best-fit shifts and shift errors in Table [3] Al- 
though many of our reference stars are from the S 32 cluster, 
comparison with the non-cluster stars does not show any large 
systematic shifts. This is as expected since typical ~ lOkm/s 
association velocities are small compared to the expected pul- 
sar space velocity, as is the velocity due to differential Galac- 
tic rotation, for the expected distances. 

Table[3]lists the best-fit frame shifts for both the raw wavde- 
tect positions and the PSF fit positions. In both cases, we com- 
pute shifts only, adopting the CXO-determined roll angle. The 
net pixel shifts are consistent between the two methods. How- 
ever, the x 2 an d frame offset errors are significantly smaller 
for the PSF-fit approach. Figure [3] shows why: the scatter of 
the PSF-fit position offsets about the best-fit line is smaller, 
with the exception of the obvious outlier Source 3. We thus 
adopt the PSF-fit measurements and offset. These finally al- 
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low us to compute the position shift (relative to the PSF-fit 
registration) of the pulsar between the two frames. The values 
for both wavdetect and improved PSF-fit detect centroiding 
are shown in Table H This table also gives the best fit radial 
(in native ACIS pixel units) and PA shifts with propagated 
errors. In Figure 4 we show the registered frames near the 
pulsar, after shifting the 2000 epoch by (69.17, 60.83) pixels 
and regridding with the CI AO dmregrid2 function. The circle 
marks the best-fit 2000 position while the ellipse marks the 
2011 epoch localization, including frame shift uncertainties. 
Between epochs the pulsar shows a significant 4.6<r shift. We 
thus have detected the proper motion of PSR Jl 809-2332; 
we comment on the implications below. 

To determine the absolute position at the 2011 epoch, we 
computed a least-squares fit using the PSF-fit positions and 
the optical positions of their USNO B-l counterparts to obtain 
a 2011.66 pulsar location of RA(2000.0)=18:09:50.25±0.03, 
DEC(2000.0)=-23:32:22.68±0.10. 



TABLE 3 
Frame Offsets: Two Fit Methods 



Method 


Coord. 


Shift (pix) 


X 2 /DoF 


Wavdetect 


X 


69.12 ±0.18 


29.9/12 


Wavdetect 


Y 


60.86 ±0.15 


16.2/12 


PSF Fit 


X 


69.17 ±0.14 


18.4/12 


PSF Fit 


Y 


60.83 ±0.09 


5.0/12 
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Fig. 3. — Residual source displacements between epochs (5 = new po- 
sition — old position), relative to the PSF-fit frame offset (Table 0. Blue 
indicates wavdetect positions and errors, while red indicates positions and 
errors obtained from the FoM maps. Source 3, an obvious outlier, is pruned 
from the frame-shift fits. The pulsar shows small but significant shifts. 



4. X-RAY SPECTRA 

For our spectral study we define extraction regions for the 
several compact and extended sources. On the smallest scales 
these are a2" point source aperture and a 3.5" x 5" 'Torus' 



TABLE 4 

Pulsar Displacement (201 1 .66 - 2000.71) 



Method X a Y a R a (radial) b 



Wavdetect -0.33 ± 0.23 -0.40 ± 0.21 0.52 ±0.15 140.6 ± 24.2 

PSF Fit -0.27 ±0.19 -0.54 ±0.16 0.60 ±0.13 153.3 ± 18.4 

a pixels (1 pixel = 0.49") 
b measured in degrees CCW from north 




Fig. 4. — Left: Chandra image of the pulsar from the 2000 observation, 
registered to the 2011 frame. Right: Matched Chandra image from the 2011 
observation. Positions and errors (black circle 2000, green ellipse 2011) are 
taken from the FoM maps. The 201 1 uncertainty includes the frame-tie errors 
added in quadrature (Table[3}. The offset is significant at ~4cr. 

ellipse, transverse to the major nebula axis. This small aper- 
ture is visible, but unlabeled on Figure 5. On this figure we 
label the larger scale spectral extraction apertures: a brightest 
nebula 'Trail', a surrounding 'Inner Nebula' and three large 
regions of low surface brightness (an 'Upstream' region ahead 
of the pulsar and Northern and Southern extensions of the 
outer nebula). In each case the apertures exclude any enclosed 
smaller scale region. As noted in Section 3 we also extracted 
point-source spectra and fit absorbed Mekal models to the ob- 
vious stellar sources to derive spectra for computing model 
PSFs. These stellar sources were not exceptional and will not 
be discussed further here. Background regions are defined on 
source-free portions of the 10, II and 13 chips; scaled back- 
grounds are subtracted from the source spectra. 

Spectra are extracted with the CI AO version 4.3 specextract 
function, which also computes response files. We group the 
spectra to a signal-to-noise/bin of 3, and fit all diffuse regions 
jointly with old and new-epoch data to an absorbed power- 
law model, using Sherpa. For the point source (pulsar) there 
is no evidence for variability, with 7.4±0.9xl0 -3 cps during 
the first epoch and 7.6 ± 0.5 x 10 -3 cps in the second. We 
therefore fit jointly to an absorbed power-law plus a thermal 
component. Table \5\ gives the fit values and 68% confidence 
errors. If we assume a common origin for the various com- 
ponents we can improve the spectral constraints by fitting for 
a global Nh- Using the brightest diffuse regions (Trail + In- 
ner Nebula) we obtain N H = 1.3 ± 0.2 x 10 22 cm~ 2 and fix 
this for the other spectral regions. The Galactic HI surveys 
(nhtool) show a column of ~ 5 x 10 21 cm -2 to d ~ 0.5 kpc 
suggesting a source distance > 1 kpc. 

For the point source (neutron star), we first attempted a joint 
BB+PL fit, but found that the PL index was driven to very soft 
values T > 6, indicating a poor fit to the thermal spectrum. If 
fixed at a more physical T = 2 the fit is unacceptable. How- 
ever, if we adopt a neutron star atmosphere model (Sherpa 
model xsnsa) we obtain a reasonable reduced x 2 an d power- 
law index. This model includes a number of parameters, three 
of which we fix: neutron star gravitation mass = 1.5M , neu- 
tron star radius =10 km, neutron star magnetic field = 10 13 
G. The remaining variables are the effective temperature and 
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Fig. 5. — Merged 0.5-7 keV image showing the spectral extraction regions 
and background regions (dotted). The torus region is too small to be clearly 
seen, but is a 3.5" x 5" ellipse centered on the pulsar and excising the inner 
2" 

model normalization. The model normalization provides an 
estimate of the pulsar distance, which for the fitted parame- 
ters yields d = 0.3^0.? kpc. This distance estimate is not very 
meaningful since, as usual, the thermal emission likely has 
contributions from a heated polar cap or soft magneto spheric 
power law. With the addition of such components the inferred 
distance increases ; it is in any ev ent consistent with the Nh 
constraint and the lOka et al.l (|1999[) distance estimate. 

The diffuse spectral indices are not atypical of hard PWN 
emission. We do not see any clear systematic softening as one 
moves farther from the pulsar. However both the 'Upstream' 
region and the 'Outer Neb -S' appear significantly softer than 
the rest of the nebula. This is especially true if Nh is allowed 
to vary, when somewhat larger T values are preferred. The 
upstream emission is difficult to explain given the detected 
proper motion of the pulsar, unless it represents a foreground 
or background structure. This leads to the tentative sugges- 
tion that the hard-spectrum PWN is superposed on a softer 
background arc composed of the 'Upstream' and 'Outer Neb 
-S' regions. This emission may be unrelated to the PWN. 
A plausible origin is a reverse shock structure in the SNR. 
Unfortunately we lack adequate S/N to confirm the spectral 
differences or even to test whether a thermal fit is more suit- 
able than a power-law model. In contrast, the 'Outer Neb -N' 
region is formally very hard, quite distinct to the Southern 
branch. Again the limited counts preclude any detailed spec- 
tral study. 

5. DISCUSSION AND CONCLUSIONS 

The displacement of 0.60 ± 0.13 pixels over 10.95 years 
gives a proper motion of /i = 27 ± 5 masyr -1 . This corre- 
sponds to a space velocity of 231 ± 46(d/1.8kpc) kins" 1 , a 
modest young pulsar space velocity. Interestingly, the proper 
motion vector points back to the birthsite of I = 7.53°, 
b = —1.68° inferred within the < y 1.5 ° diameter radio shell 
G7.5-1.7 bv lRoberts & Broganl (120081) . This is 1300" from 



the pulsar so, taken at face value, our proper motion implies 
a kinematic age of 48 ± lOkyr. Given the characteristic age 
r c = 68 kyr, one gets an initial spin period Po ~ 104 ± 20 ms, 
for a braking i ndex n = 3. 

As noted by Rob erts & Bro gan ( 2008), the pulsar projected 
offset is only ~ 45% of the SNR radius, so if the PWN lies in 
the SNR interior it may be just passing from being 'crushed' 
by the reverse shock to forming a w ell-defined bow shock 
(Ivan der Swaluw. Downes & Keegan 2004). Accordingly, we 
should not b e too surprised that the morphology at the PWN 
apex (Figure |l(b)| ) is unclear. Certainly the bulk of the hard 
spectrum emission trails the pulsar, following the proper mo- 
tion axis and the direction to the explosion site. However, 
there are significant counts bracketing the point source and 
extending ~ 5" transverse to the pulsar motion. Since the 
pulsar spin axis appe a rs to correlate with t he proper motion 
(iJohnston etaLl 120051: iNg & RomanH 120071) such transverse 
extension tends to be equatorial. It is tempting to infer that the 
blocky PWN head is the result of an anisotropic pulsar wind, 
concentrated in an equatorial torus, with a spin axis/Earth 
line-of- sight angle ( > 75°, i.e. a torus viewed at large in- 
clination angle. It is interesting to compare with predictions 
for the observed 7-ray pulse, which is a fairly narrow dou- 
bl e with peak separat i on A = 0.35. Examining the 'Atlas' 
of Romani & Watters (2010), we see that 7-ray pulsars with 
A « 0.35 should have 70° > ( > 80° and magnetic axis in- 
clination a < 60° for outer-magneto sphere dominated emis- 
sion. There is little phase space for such pulsars to be radio de- 
tected. TPC-type models produce such A for a wide range of 
C < 70°, but many of the solutions should be radio detected. 
Thus both model classes are allowed, although an outer-gap 
type interpretation seems preferred. A good measurement of 
£ from a detailed map of the PWN head could check this in- 
terpretation, but would require a rather long ACIS exposure. 

Looking ahead of this transverse structure, we see that the 
hard spectrum diffuse emission has a fairly abrupt cutoff at 
~ 1" (Figure lb). We can compare this with the standoff 
angle for an isotropic ram-pressure confined pulsar wind 

°° 8 - (i^) ' - 1 - 3 " n ~ 1/2 d - (1) 

for our measured proper motion, an ambient number density 
n cm -3 and a distance 1.8 kpc. It seems that a small standoff 
is not unexpected, but given the apparent anisotropy in the 
wind momentum and likely anisotropy in the SNR interior, 
no strong conclusions should be drawn. 

Of course there is emission 'Upstream', ahead of the pulsar 
motion. The chip gap spanned this region in the 2000 data 
so it was difficult to draw morphological conclusions. In the 
combined data it seems that this emission is morphologically 
distinct to, and substantially softer than, the PWN trail. This 
larger-scale soft emission is also seen in archival XMM and 
ASCA data where it also appears distinct from the harder PWN 
trail. In our data it appears to connect to the "Outer Neb - S" 
region, so we can posit that this extended softer emission may 
be a background or foreground structure associated with the 
host SNR. If the absorption is left free, these regions seem to 
prefer a slightly higher Nh than the central PWN. Given the 
extensive patchy molecular gas and obscuration in this region, 
these small differences are not particularly telling. 

We have examined a new Chandra ACIS exposure of PSR 
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TABLE 5 
X-Ray Spectral Fits 



Region 




r/LogT eff 


Abs. Flux b 


Unabs. Flux b 


x 2 /d.o.f. 


Trail+Inner 


qi+0.20 

1 - ol -0.18 


1.83 ±0.16 


6.3±i;i 


n.9±i:I 


112.9/199 


PSR 


1.3 C 


2.00 ±0.54/ 5.9 ±0.1 






12.7/26 


Torus 
Trail 

Inner Neb 
Upstream 
Upstream 
Outer Neb 
Outer Neb - S 
Outer Neb - S 
Outer Neb - N 


1.3 C 
1.3 C 
1.3 C 
1.3 C 

1 8 + - 5 

1.3 C 
1.3 C 
1.8 C 
1.3 C 


--+0.41 
L - ' ' -0.38 
-, oo + 0.19 
J--OO-0.18 

1.79 ±0.09 

9 91 +0.19 
Z - Z1 -0.18 

Z - DD -0.41 

1.56 ±0.10 
2.19 ±0.14 
2.61 ±0.16 
1.00 ±0.17 


n q+0.2 

°- 3 -o.i 
0.9 ±0.2 

2.3 ±0.4 

o 9 + 1.8 
z - z -0.9 
o 4 +1.0 
8 - 4 -0.9 

3.9 ±0.5 

o 7 +0.6 
6( -0.5 

4 '^-0.9 


n «+0.3 
°- 6 -0.2 
1 8 + - 4 

«- 4 -0.8 

5 4+1-0 

D - 4 -0.9 

12.9^° 4 4 

9.0+1? 
13.3 + ?-g 
5.811:1 


5.3/8 
14.9/29 
76.3/159 
93.0/152 
91.6/151 
196.6/348 
108.2/209 
113.0/209 
93.5/178 



interstellar absorption x 10 22 cm 2 
6 0.5 — 7keV fluxes in units of 10 -13 erg cm -2 s _1 
c held fixed 



Jl 809— 2332 to study the fine scale X-ray morphology. Com- 
parison with the 2000 exposure yields a > 4<j detection of 
the pulsar proper motion and supports the association with 
SNR G7.5-1.7. The PWN has slightly extended emission at 
the apex, somewhat larger than expected from a bow shock. 
While this suggests that the pulsar wind is concentrated trans- 
verse to its velocity, the possibility of anisotropies in the SNR 
interior discourages strong conclusions and such disturbed 
morphology is not unexpected. Overall, however, these new 
data support a basic picture for the PSR Jl 809— 2332 system: 
a ~ 50ky pulsar viewed near the spin equator, but with small 



magnetic inclination so that the radio beam misses the Earth. 
The pulsar, traveling at ~ 230km s _1 , is followed by a trail of 
hard PWN X-ray emission and is approaching the outer region 
of its composite host SNR located at d < 2 kpc. 
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